###Exercicios - Aula 6

getwd()
setwd("C:/Users/Edmundo/Documents/Meus Documentos/BIANCA/MESTRADO/Linguagem R")
getwd()

##Crie seus dados

help(rnorm)
conjA <- rnorm(10,mean=6,sd=3)
conjA
conjB <- rnorm(10,mean=7.5,sd=3.2)
conjB
#
source("simula.r")
dif.conjs <- simula(conjA,conjB,nsim=2000)
dif = abs(mean(conjA)-mean(conjB))
prob <- sum(dif.conjs>=dif)/length(dif.conjs)
prob
#P=0.848. A hiptese nula de que A=B no pode ser rejeitada.
dif.conjs2 <- simula(conjA,conjB,nsim=2000,teste="uni")
dif2 <- abs(mean(conjA)-mean(conjB))
prob2 <- sum(dif.conjs2>=dif2)/length(dif.conjs2)
prob2
#P=0.4205. Tambm no podemos rejeitar a hiptese de que A > B.
t.test(conjA,conjB)
P=0.8314

##Caixetas de novo?!

caixeta <- read.table("caixeta.csv",as.is=TRUE,head=TRUE,sep=",")
caixeta
#
caixeta$local=factor(caixeta$local)
caixeta$parcela=factor(caixeta$parcela)
caixeta$fuste=factor(caixeta$fuste)
caixeta$arvore=factor(caixeta$arvore)
#
raio <- caixeta$cap/2
raio
caixeta$area.basal <- pi*(raio^2)
caixeta$area.basal
#
help(aggregate)
ABamostra <- aggregate(caixeta$area.basal,list(caixeta$local,caixeta$parcela),FUN=mean)
ABamostra
names(ABamostra) <- c("parcela","local","area")
#
boxplot(caixeta$area.basal~caixeta$fuste)
boxplot(ABamostra$area~ABamostra$local)
#Valores da tabela de ANOVA:
#FUSTES#
media.t <- mean(caixeta$area.basal)
#Desvio quadrtico total:
ss.tot <- sum((caixeta$area.basal-media.t)^2)
ss.tot
#Intra-grupos:
m.chauas <- mean(caixeta$area.basal[caixeta$loca=="chauas"])
m.jureia <- mean(caixeta$area.basal[caixeta$loca=="jureia"])
m.retiro <- mean(caixeta$area.basal[caixeta$loca=="retiro"])
ss.chauas <- sum((caixeta$area.basal[caixeta$local=="chauas"]-m.chauas)^2)
ss.jureia <- sum((caixeta$area.basal[caixeta$local=="jureia"]-m.jureia)^2)
ss.retiro <- sum((caixeta$area.basal[caixeta$local=="retiro"]-m.retiro)^2)
ss.intra <- ss.chauas+ss.retiro+ss.jureia
ss.intra
#Entre grupos:
ss.entre = ss.tot - ss.intra
ss.entre
#Graus de liberdade:
GLtot <- length(caixeta$area.basal)-1
GLtot
GLintra <- length(levels(caixeta$local))-1
GLintra
GLentre = GLtot - GLintra
GLentre
#Desvios mdios:
DMentre = ss.entre/GLentre
DMintra = ss.intra/GLintra
#
F = DMentre/DMintra
F
pf(F,GLentre,GLintra,lower.tail=FALSE)
#PARCELA#
media.parcela <- mean(ABamostra$area)
#Desvio quadrtico total:
ss.tot.PAR <- sum((ABamostra$area-media.parcela)^2)
ss.tot.PAR
#Intra-grupos:
ss.chauas.PAR <- sum((ABamostra$area[ABamostra$local=="chauas"]-mean(ABamostra$area[ABamostra$local=="chauas"]))^2)
ss.jureia.PAR <- sum((ABamostra$area[ABamostra$local=="jureia"]-mean(ABamostra$area[ABamostra$local=="jureia"]))^2)
ss.retiro.PAR <- sum((ABamostra$area[ABamostra$local=="retiro"]-mean(ABamostra$area[ABamostra$local=="retiro"]))^2)
ss.intra.PAR <- ss.chauas.PAR+ss.retiro.PAR+ss.jureia.PAR
ss.intra.PAR
#Entre grupos:
ss.entre.PAR = ss.tot.PAR - ss.intra.PAR
ss.entre.PAR
#Graus de liberdade:
GLtot.PAR <- length(ABamostra$area)-1
GLtot.PAR
GLintra.PAR <- length(levels(ABamostra$local))-1
GLintra.PAR
GLentre.PAR = GLtot.PAR - GLintra.PAR
GLentre.PAR
#Desvios mdios:
DMentre.PAR = ss.entre.PAR/GLentre.PAR
DMintra.PAR = ss.intra.PAR/GLintra.PAR
#
F.PAR = DMentre.PAR/DMintra.PAR
F.PAR
pf(F.PAR,GLentre.PAR,GLintra.PAR,lower.tail=FALSE)
#Funo ANOVA:
anova.fuste <- aov(area.basal~local,data=caixeta)
summary(anova.fuste)
anova.parcela <- aov(area~local,data=ABamostra)
summary(anova.parcela)
#
ss.entre/ss.tot
#Local explica 6.43% da variao
ss.entre.PAR/ss.tot.PAR
#Local explica 100% da variao (???)